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1. Introduction. Statement of the Problem 

In order to find an approximate representation of a function ip{x) by elements of a certain finite 
dimensional space, it is possible to use values of this function at some finite set of points Xk, k = 
1,2,...,N. The corresponding problem is called the interpolation problem, and the points the 
interpolation nodes. 

There are polynomial and spline interpolations. Now the theory of spline interpolation is fast 
developing. Many books are devoted to the theory of splines, for example, Ahlberg et al [1], de Boor 
[2], Schumaker [3], Laurent [1], Attea [5], Stechkin and Subbotin [6J, Vasilenko [7|, Arcangeli et al 
[8], Ignatov and Pevniy [9], Korneichuk et al [lOj, in the numerical analysis literature or Wahba jll], 
Eubank [12j . Green and Silverman [13j, Berlinet and Thomas- Agnan [14J in the statistical one. 

If the exact values (p{xk) of an unknown smooth function (p{x) at the set of points {xk, k = 
1,2, A^} in an interval [a,b] are known, it is usual to approximate ip by minimizing 



(1.1) / {g^"'Hx))'dx 



in the set of interpolating functions (i.e. g{xk) = ^{xk), k = 1,2, A^) of the Sobolev space L^\a, h). 
It turns out that the solution is a natural polynomial spline of order 2m with knots xi, ...,xm called 
the interpolating D™' spline for the points {xk,ip{x]S)). In non periodic case first this problem was 
investigated by Holladay [15] for m = 2 and the result of Holladay was generalized by de Boor [16] for 

any m. In the Sobolev space l'^'' of periodic functions the minimization problem of integrals of type 
p.lj) were investigated by I.J.Schoenberg [17] and M.Golomb [18] . 

In the present paper we deal with optimal interpolation formulas. Now we give the statement of 
the problem of optimal interpolation formulas following by S.L. Sobolev [19]. It should be noted that 
connection between the minimization problem of integrals of type p.ip and the problem of optimal 
formulas was shown in [20J. 

1 



2 



KH.M. SHADIMETOV, A.R. HAYOTOV, N.H. MAMATOVA 



First we recall the definition of Sobolev space (0,1) of periodic functions (see HUES]). 
Suppose in the set M of real numbers the function ip{x) has local summable derivatives up to order 

m, and also for the interval [0, 1] the integral {^p^'^\x)) dx is bounded. Assume that the function 
'^{x) is periodic: (p{x + 7) = '^{x), x G M, 7 G where Z is the set of integer numbers. 

Every element of the space 4"^ (0,1) is a class of functions which differ from each other by constant 

term. The norm of functions in the space L^\{), 1) is defined by formula 

1/2 



ip'->{x)^ 

Now following [19j we consider interpolation formula of the form 




dx 



[1.2) 



(fix) 



N 

Y. 

k=l 



Ck{x) ■ ip{xk) 



in the space l'^\o, 1). Here points Xk G [0, 1] and parameters Cfc(x) are respectively called the nodes 
and the coefficients of the interpolation formula ()1.2p . 

The difference (p — P^p is called the error of the interpolation formula (|1.2p . The value of this error 
at some point z is the linear functional on functions (p, i.e. 



N 



?, if) = ip{z) - P^{z) = ip{z) - J2 Ck{z)ip{xk) 



k=l 



(1.3) 




TV 



6{x - z) - ^ Ck{z)6{x -Xk)\ * (/>o(x) 



k=l 



ip{x)dx, 



where 6{x) is the Dirac delta-function, (po{x) = '^S{x — /3), here /? takes all integer values and 

/3 



:i.4) 



e= (^6{x-z)-Y^Ck{z)5{x-Xk)^ 



is the error functional of interpolation formula (|1.2p and belongs to the space 4"^* (0,1). The space 

4"^* (0,1) is the conjugate space to the space -L2™'*(0, 1) and consists of all periodic functionals ([I 
which are orthogonal to unity, i.e. 

(1.5) (^,1) = 0. 

By Cauchy-Schwartz inequality 

i(£,v.)i<iiv.i4"^(o,i)ii-ii^i4'"^*(o,i)ii 
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the error p.3p of formula ()1.2p is estimated with the help of the norm 



£i4™^*(o,i) 



sup 



of the error functional (II. 4p . Consequently, estimation of the error of interpolation formula (ll.2p 

on functions of the space l'^^ (0, 1) is reduced to finding the norm of the error functional £ in the 

conjugate space l'^^*{0, 1). 

Therefore from here we get the first problem. 

Problem 1. Find the norm of the error functional £ of interpolation formula II LS^) in the space 



L 



(m)* 



(0,1). 



Obviously the norm of the error functional £{x) depends on the coefficients Ck{z) and the nodes 
Xk- By optimal interpolation formula we call such formula which the error functional in given number 

of the nodes has the minimum norm in the space L2^^*[0, 1). If the nodes Xk are the points of a 
lattice, i.e. are located at points of the form Xk = hk then such interpolation formula is called the 
lattice interpolation formula. Here h is small parameter and is called the step of the lattice. 

The main goal of the present paper is to construct the lattice optimal interpolation formula in the 

space l'^\o,1) for the nodes Xk = hk, i.e., to find the coefficients Ck{z) satisfying the following 
equality 



(1.6) 



^(x)|4-)*(0,l) 



inf 



£{x)\Lt^*{0,l] 



Thus in order to construct the lattice optimal interpolation formula in the space l'^\o, 1) we need 
to solve the next problem. 

Problem 2. Find the coefficients Ck{z) which satisfy equality when the nodes are located in the 
lattice, i.e., Xk = hk. 

In the present paper the lattice optimal interpolation formulas are constructed in the Sobolev space 

L^^\o,l) of periodic functions. First such problem was stated and investigated by S.L. Sobolev in 

|19j . where the extremal function of the interpolation formula was found in the space W2"^\ 

This paper is organized as follows. In Section 2 we determine the extremal function which cor- 
responds to the error functional £{x) and give a representation of the norm of the error functional 
p.4j) . Section 3 is devoted to minimization of ||^|| with respect to the coefficients Ck{z). We ob- 
tain a system of linear equations for the coefficients of the optimal interpolation formula in the space 

L^\o,l). Moreover, the existence and uniqueness of the corresponding solution is proved. In Sec- 
tion 4 we consider the problem of construction of lattice optimal interpolation formulas. In Section 5 
we prove some new properties of the discrete analogue of the differential operator d^^/dx^"*. Using 
these properties explicit formulas for coefficients of lattice optimal interpolation formulas are found 
in Section 6. In Section 7 the norm of the error functional of lattice optimal interpolation formulas is 
calculated. Integrating obtained lattice optimal interpolation formula the optimal quadrature formula 
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in the space L'i^\o, 1) is obtained in Section 8. FinaUy, in Section 9 we give formulas for coefficients 
which are very useful in practice and we present some numerical results. 

2. The extremal function and representation 
of the norm of the error functional 

In this section we solve Problem 1, i.e. we find explicit form of the norm of £{x). 

For finding the explicit form of the norm of the error functional i{x) in the space 4™) (0,1) we 
use concept of its extremal function which was introduced by S.L.Sobolev [191 [21]. The function u{x) 

from L2^\o, 1) is called the extremal function for the error functional £(x) if the following equality is 
fulfihed 



4")* (0,1) 



u 



4-) (0,1) 



The space L^"^(0,1) is Hilbert space and the inner product in this space is given by formula 



X. 



According to the Riesz theorem any linear continuous functional i{x) in a Hilbert space is represented 
in the form of a inner product 

(2.1) i£,ip) = {iPe,ip)m 

for arbitrary function ip{x) from L^\o, 1). Here ipi{x) is a function from l'^\o, 1) is defined uniquely 
by functional l{x) and is the extremal function. Integrating by parts the expression in the right hand 
side of (j2.ip and using periodicity of functions ip{x) and Tpi{x) we get the following equality 



(A 9^) 



'2m 



'4)t{x) ■ ip{x)dx. 



Thus the function ipi{x) is the generalized solution of the equation 



(2.2) 

with the boundary conditions 



a 



0,2m - 1. 



For the extremal function the following holds 

Theorem 2.1. Explicit expression for the extremal function ipe{x) of the error functional ji.^l ) is 
defined by formula 

N 

(2.3) ^^(a;) = (-l)- B2m{x-z)-Y,Ck{z)-B2Ux-Xk) + do 

k=l 

where B2m{x) = Yl, ^^{2ni^)^"^^^ ^'^ Bernoulli polynomial, do is a constant. 
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Proof. Here we use the following formulas of Fourier transformations given in [21] 
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F[ip{x)] = / ex.p{2iTipx) ■ (p{x)dx, F ^[^{p)] = / ex.p{—2-7Tixp) ■ (p{p)dp 



Convolution of two functions is defined by formula 



oo 

f{x)*g{x)= J f{x-y)-g{y)dy. 



Applying to both sides of equation (|2.2|) the Fourier transformation and known formulas 
F[6{x - z)] = e2^^P^ F[(I)q{x)] = Mp) (see [21]) we get 



(2.4) {27Tipy^F[Mx)] = (-1) 



N 



exp(27ripz) - ^ Cfc(z) exp(27ripxfc) </>o(p) 



k=l 



By virtue of (jl.Sp the right hand side of ()2.4p is zero at the origin. Therefore we can divide both 
sides of ()2.4p by (27rip)^™. The function F[ip£{x)] is defined from equation (|2.4p up to the following 
expression 



2m.- 1 



(-l)"^do<5(p) + d^D-6{p). 



a=l 



But as known a periodic solution of the homogenous equation corresponding to equation ()2.2p is a 
constant term then all terms except (— l)™(io(5(p) should be omitted. Thus from (j2.4p we get 



(2.5) F[Mx)] = {-irdo5{p) + 



. V , . X E Cfc(2;)exp(27ripj;fc) (?:)o(p) 
exp(27rzp2;) (/>o(p) fc=i 



(2vrp) 



2m 



(2^p) 



2m 



Changing the function (f)Q{p) by series of 6 functions and applying the inverse Fourier transformation 
to both sides of (|2.5p we get (|2.3p . Theorem 12.11 is proved. □ 



Now we obtain representation for the norm of the error functional i. Since the space l'^''\o, 1) is 
the Hilbert space then by the Riesz theorem we have 



(2.6) 



i 



4™)* (0,1) 



4™) (0,1) 



4™)* (0,1) 
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Using formulas ([Ll]), (I23]), <^Ml we get 



e{x) 



4™^* (0,1) 



N 



X yB2m{x - Z 

Hence, using the condition (ll.5p we have 



i{x)M^)dx ={-ir J \^5{x - ^) - E ^'^(^) • '^(^ - ^fc) j * Mx) 

N \ 
Ck{z) ■ B2m{x - Xk) + do 
k=l J 

N \ 

'o(^) ~ XI Ck{z)5{x - Xk) * (poix) j 

A;=l / 
N \ 

Ckiz) ■ B2U^ - Xk) dx + (-l)"(io {l{x), 1) . 

k=l J 



X B2ra{x - Z 



(-l)™y \5{X-Z)* 






4")* (0,1) 


2 1 


oo 


l{x) 


= (-ir/ 


^ / 5{x - z - y)6{y - f3)dy 









. P - oo 



/ ^(^ - - y)Siy - /3)dy 
k=l -'^ 



N 



B2m{x - z) - Ck{z) ■ B2m{x - Xk) 



k=l 



dx. 



Taking into account definition of Dirac's delta-function we obtain 



i{x) 



4"^)* (0,1) 



N 



Y,5{x-z-p)-Y, Ck{z)5{x -Xk-P) 



k=l 



i-iy 



X ^B2m{x - z) -"Y^^ Ck{z) • (x " 2;^)^ dx = 

1 TV 1 

/ Y^i^ ~ ^ ~ (^) ' B2m{x - z)dx- ^ Ck{z) / X 5(x - z - /3) • B2m{x - Xk)dx- 

0/3 k=i i p 

} 

-'^Ckiz) / y^d{x - Xk - f3) ■ B2m{x - z)dx+ 

k=i i 8 



N N 



+ 



Cfc(2:)Cfc'(z) X / 5(X - Xfc - f3)B2m{x - Xk')dx 
k=lk'=l P Q 
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Hence with the help of the characteristic function X[o,i](^) of the interval [0, 1] square of the norm of 
the error functional (jl.4p we reduce to the form 



iix) 



4™)* (0,1) 



(-1) 



P - CO 

oo 

^Cfc(z)^ / X[Q,i\{x)5{x - Z - f3) B2m{x - Xk)dx 
k=l P 

N °9 

-^C'fc(z)^ / X[o,i]{x)d{x - Xk - l3) B2m{x - z)dx 

N N °9 
+ ^ ^ Cfc(2)Cfc/(z)^ / X[0,l]ix)Kx - Xk - /3) B2m{x - Xk)dx 
k=l k'=l B _^ 



N 



From here applying definition of Dirac's delta-function we have 



l{x) 



4™)* (0,1) 



(-ir Emi](^+/5)^2m(/3) 

k=l (5 
N 

X Ck{z) X X[0,1] {xk + P)B2rn{xk + P - z] 



k=l 
N 



N 



+ X] ^''^^^ X] ^k'{z)^X[0,l](.Xk + l3)B2m{Xk + P -Xk' 
k=l k'=l /3 

Taking into account that Xk G [0,1], z £ [0,1], X^/j X[o,i] (?/ + Z^) = 1 (see. [21J), y G [0,1], and 1- 
periodicity, symmetry of B2m{x), that is B2m{x + 7) = B2mix), B2m{—x) = B2m{x) from the last 
equality we get the following representation of ||^(a;)||^ 



l{x) 



4"^)* (0,1) 



(-1)' 

N 



N 



(2.7) 



B2ra{^)-2Y,Ck{z)B2ra{z-Xk) + 
k=l 

N 

+ X X '^k'{z)B2m{Xk - Xk') 



k=l k'=l 



Thus Problem 1 is solved. 

Further in next sections we solve Problem 2. 
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3. Existence and uniqueness of optimal interpolation formula 

For finding the minimum of the norm (j2.7p under the condition (jl.Sp we use the Lagrange methods. 
For this we consider the function 

^{C{z),X) = \\if + 2{-irX {e{x),l), 

where C(z) = {Ci{z),C2{z), ...,Cn{z)). 

Equating the partial derivatives of 'if{C{z), A) by Ck{z), {k = 1, 2, N) and A to zero we get the 
fohowing system of hnear equations 

N 

(3.1) ^^Ck{z)B2m{Xk' - Xk) + \ = B2m{z - Xk'), /c' = 1, 2, iV, 
k=l 

N 

(3.2) Y.'^k{z) = l. 

k=l 

System (3.1)-(3.2) has a unique solution and this solution gives the minimum to ||^(x)|p under the 
condition ()3.2p . 

It should be noted that uniqueness of the solution of such type systems were also investigated in 

[ZlElEllZIllSlj. 

The uniqueness of the solution of system (I3.ip - (l3.2p is proved following [221 Chapter I]. For com- 
pleteness we give it here. 

First in (12. 7|) we change of variables Ck{z) = Ck{z) + dk{z) then (12. 7|) and system (I3.ip - (l3.2p have 
the following form 



(-ir 



N N N 



'^Ck{z) ^ Ck'{z)B2m{xk - Xk') - 2^(Cfc(z) + dk{z))B2m{z - Xk) 



_k=l k'=l k=l 

N N 

(3.3) (.dk{z)dk'iz) + Ck{z)dk'{z) + Ck'{z)dkiz))B2miXk - Xk') + B2m 

k=l fc'=l 

N N 

(3.4) ^ Ck{z)B2m{Xk' -Xk) + \ = B2m{z - Xk') - ^ dk{z)B2m{xk' " Xk), 
k=l k=l 

k' = 1,2,...,N, 

N 

(3.5) Y,Ck{z) = 0, 

k=l 

where dk{z) is a partial solution of the equation (j3.2p . 

Hence we directly get that the minimization of (2.7) under the condition (j3.2p with respect to 
Ck{z) is equivalent to the minimization of expression ()3.3p with respect to Ck{z) under the condition 
(|3.5p . Therefore it is sufficient to prove that system (j3.4p - (j3.5p has a unique solution with respect 
to unknowns C{z) = {Ci{z),C2iz), ...,C m{z)), A and this solution gives the conditional minimum for 

From the theory of conditional extremum it is known the sufficient condition in which the solution of 
system (I3.4p - (j3.5p gives the conditional minimum for ||^||^ on manifold ()3.5p . It consists of positiveness 



OPTIMAL INTERPOLATION FORMULAS 

of the fohowing quadratic form 



N N 



(3-6) ^(C(.)) = E E 7^7W7\^^(^)^^'(^) 

t^i^i dCk(z)dCk'{z) 

on the set of vectors C{z) = {Ci{z),C2{z), ...,C n{z)) under the condition 

(3.7) SC = 0, 

where S = (1, 1, 1) is the N component vector. 

In our case this condition is fulfilled, i.e. the following holds 

Lemma 3.1. For any non zero vector C{z) = {Ci{z),C2{z), ...,C ]\f{z)) G lying in the subspace 
|g. 7|) the function ^{C{z)) is strictly positive. 



Proof. From definition of the function ^'(C(z), A) and (13. 6p it follows that 

N N 

(3.8) <!>(C{z)) = 2Y,^(.-^rB2m{xk-Xk')Ck{z)Ck'{z). 

k=i k'=i 

We consider the functional 

N 



%(2)(^) = V2'^Ck{z)S{x - Xk) * (t>o{x). 



k=l 

By virtue of condition (j3.7p the functional ^c(2)(^) belongs to the space L2™'**(0, 1). Thus this func- 
tional has the extremal function 11^(2) ^ ^2"^^^ ^) which is the solution of the equation 

(3-9) ^-C(.)(-) = (-ir%(.)(-)- 

As 

ttc(z) (^) following function 

N 
k=l 



Square of the norm of 'Uc(z)(^) coincide with the function ^{C{z)) in the space 4^^(0,1), i.e. 

, N N 

(3.10) \\u-^^^^ix)\L^r\0A)\\' = (^ciz)(^^ 



k=i k'=i 



From here taking into account ()3.8p we conclude that for non zero C(z) the function $(C(z)) is strictly 
positive. Lemma |3. II is proved. □ 



If the system (j3.4p - (j3.5p has a unique solution then the system (j3.ip - (|3.2p has also a unique solution. 
The following holds 

Lemma 3.2. The main matrix Q of the system ^3.4\ )- f375\} is non singular. 
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Proof. The homogenous system corresponding to system ()3.4p ~ ()3.5p have the following matrix form 



(3.11) 



Q 



C{z) 
X 



where B is the matrix with elements h 



13 



B 1 
S 

B2miXi-Xj), i 



C{z) 
X 



0, 



1, 2, N, j = 1, 2, N, S is the vector 



Ix N and S = (1,1,..., 1). 

We show that system (j3.1ip has the unique solution C(z) = (0,0, ...,0), A = 0. 

Assume C{z), X is a solution of ()3.1ip . Since a solution of equation (|3.9p is determined up to 
constant term then as Uqi^z)^''^) take the following function 

^C(.)(^) = (-l)"\/2 Cfc,(z)52™(x -Xk>) + X\. 

\k'=i J 

But from (j3.1ip it is clear that Uc(z)i-^k) = 0. Then for the norm of the functional ^c(z){x) we have 



4'")*(o,i) 



N 



on the other hand from (j3.10p we get 



^c(^)(^)'^c(^)(^)) =5^Cfc(z)-nc(^)(xfe) = 



k=l 



Lf)*(0,l) 



N N 



%(^)(x),Uc(^)(x)j = 2j2Yi-'^rB2m{Xk - Xk')Ck{z)Ck'{z) 



k=l k'=l 



which is possible only when Ck{z) = 0, k = 1,N. 

Hence by virtue of (j3.1ip we obtain A = 0. Lemma 13.21 is proved. 



□ 



From ()2.7p and Lemmas 3.1, 3.2 it follows that in fixed values of the nodes Xk square of the norm 
of the error functional £{x) being quadratic function of the coefficients Ck{z) has a unique minimum 



in some concrete value Ck{z) = Ck{z). Interpolation formulas with the coefficients Ck{z) {k = 1,N), 
corresponding to this minimum in fixed values of the nodes Xk is called the optimal interpolation 



formulas and Ck{z) {k = 1, A^) are called the optimal coefficients. 



4. The system for the coefficients of lattice optimal interpolation formula 

We consider system ()3.ip - ()3.2p from Section 3 on one dimensional lattice, i.e., suppose that the 
nodes x^^ = /17, where /i = 7 = 1, 2, . . . , A^. Below for convenience we denote [7] = /17. Then 
such an interpolation formula we call the lattice interpolation formula. Moreover in this case system 
p.ip - (j3.2p takes the following form 

N 

(4.1) Y.Cib],^)B2^[P-l\ + X = B2rn{z-m, /3 = 1, 2, iV, 

7=1 
N 

(4.2) 5](7([7],^) = L 

7=1 

Later we use the theory of discrete argument functions. The theory of discrete argument functions 
was investigated in |21j . 
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The convolution of two discrete argument functions is given by formula (cf. ^21j ) 

oo 

m*m= E /w-5[/3-7]. 

7=— oo 

Below for convenience instead of the sum ^^^^=-00 write 

Using the discrete characteristic function X[o,i][/5] of interval [0,1] and taking into account 
the definition of convolution of two discrete argument functions system (|4.ip - (j4.2p we rewrite in the 
following convolution form 

(4.3) B2m[fi]^[cm.z)-Xm[P\)+^ = B2m{z-m, [/3]G(0,1], 

N 

(4.4) E(7([/3],z) = l, 



where [p] = hp = ^. 

Now we have the following problem. 

o 

Problem 3. Find a discrete function C{[P\tz) and unknown constant A which satisfy the system 

In the solution of Problem 3 the main role plays some new property of the discrete analogue D^^^ [/3] 
of the differential operator d?^ /dx^"^. It should be noted that the properties of the discrete analogue 
of the polyharmonic operator A™ were investigated by S.L. Sobolev [21. .25j. But here we need some 
new properties of the discrete operator -D^™^ [/3] • The next section is devoted to investigation of these 
properties. 

5. Some new properties of the discrete operator L'[™^[/3] 

In the work [24j was constructed the discrete analogue Z?^*" ■*[/?] of the differential operator j^^^ . 
The discrete operator 1?^™'' is the solution of the following difference equation 

(5.1) Df)[/3]*G™[/3]=J[/3]. 



11 o\Am— L 

Here Gm[/3] = 2-{2m-i)\ ' '^I'^] equal to 1 when [j3] =0 and is equal to zero when [j3] / 0. The discrete 
operator -D[™'^ which satisfies equality (|5.ip has the following form 



( m-1 ,-. -,2m + l \fi\ 

qk-E2m-l{qk) ' - ^' 

fe = l 

fc=i 

K = l 



where < 1 are the roots of the Euler-Frobenius polynomial E2m-2i^) of degree 2m — 2. 
The definition of the Euler-Probenius polynomial is given, for example, in |22j . 
The following holds 
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Lemma 5.1. The function D^\(3] is the solution of the equation 

(5.3) hDi^\^]*B2mm = m-h, 

where B2m[P] = E ^^^felj^^, 

77^0 

(5.4) $[/3] = j;<^[/3-7/i-'], 

7 

here 5[fi — /i^^7] is equal to 1 when /? = jh"^ and is equal to zero when /3 ^ jh"^ . 



Proof. We investigate the solution of equation (j5.3p . For this we use the Fourier transformation. For 
convenience first we go from functions of discrete argument to harrow-shaped functions. According to 
the definition of harrow-shaped functions (see j2HI22]) for arbitrary function ^[13] of discrete argument 
its harrow-shaped function has the form 

^(x) = h^^J[/3]6{x - hf3). 

For <^>[/3] which defined by formula ()5.4p we have 

Y(x) = ^/i$[/3](5(x-/i/3) = ^/i^(5[/3-/i-S]^(2;-/i/3) 

= ^M(x - 7) = /i^o(a;), 
7 

where (poix) = X] "^(^ ~ 7)- 

7 

Taking into account the last equality and the identity 6{hx) = h'^^d^x) for equation (|5.3p in the class 
of harrow-shaped functions we have 

(5.5) h Dh (x)* B 2rn{x) = h(po{x) - h(po{xh ). 

The Fourier transformation of the functions (poix) and (j)oixh~^) are respectively given by formulas 



00 

(5.6) F[cf)o{x)]= y"e2^^P^-^5(x-/3)dx=^e 



2-iTipf3 

/3 



00 

(5.7) F[Mxh-^)]= I e^-^P^ J];<5(x/i-i -/3)dx =/i^e2" 

Hence using the following known formula (see [21]) 

(5.8) /i^e^"*^'^^ = ^5(p-/i-V) 

equalities ()5.6p . ()5.7p we reduce to the form 

(5.9) i^[0o(x)] = </.o(p), 

(5.10) F[(t)Q{xh-^)] = ^5(p- /i"V)c^a; =/i<^o(M- 



iphfj 
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Now we calculate the Fourier transformation of the function 

^ ( \ ( 1 ^m ^ exp(-27ri7/t/3)^(x - 
B2^{x) = {-1) 2^/1 2^ f 271-7) 2™ • 

By definition of the Fourier transformation we have 

00 

pr^ ^T^l - [ p2-ipx ^ r-rW-r - \^ /, exp(-27ri7fe/3)exp(27rip/i/3) 

^ exp(27ri/i/?(j9 - 7)) 

By virtue of equality (jS.Sp we obtain 

/i ^ exp(27ri/i/3(p - 7)) = ^ 5(p - 7 - /i" V)- 

Therefore 

F[ B 2m{x)] = J—^ 2^ 2^ -2^;^ . 

Hence, setting 7 + /i"-"^/? = /c, j = k — h~^f5, we get 

Now applying the Fourier transformation to both sides of equation (jS.Sp . using (|5.9p . (|5.10p . (|5.1ip 
and obtained equality dividing by h we get the following equation 

It is known that 



^ (7-/,-l/3)2m ^ (p_/,-l^)2„^• 
7h^Z 



By virtue of equality (jS.lSp equation (j5.12p is equivalent to the following equation 
(5.14) F\D^ W](r^^(p)) =1, v + hr^i. 



where 



(5.15) r("^)(p) 



(27r)2™ ^ /i-l^)2m 

The function rjj™'^(p) is periodic with respect to p with period /i^^, real and analytic in all p 7^ 

From equality (I5.14p we have 
(5.16) (x)] = ri:")(p). 
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Applying to both sides of equality (15.160 the inverse Fourier transformation and going from harrow- 
shaped functions to functions of a discrete argument we find (|5.2p (see [24j ) . This completes the proof 
of Lemma 15. 11 □ 



Lemma 5.2. For the convolution of discrete functions D^^\l3] and ex.p{27riah/3) the following holds 



(5.17) Dj^""^ [/?] * exp(2^ia/i/3) = (-l)'^(27r)2" exp(2^i(j/i/3) 



-I -1 



E 



(7 - cr/l)2"^ 



, ah Z, 



(5.18) 



Mr 



exp(27ricr/i/3) =0, ah e Z. 



where Dj^™^[/3] is defined by equality \5. 



Proof. The convolution of two harrow-shaped functions D h (^) exp(27ric7x) we denote by T (a;), 
i.e 

T {x) = D h i^)* exp{2iTiax). 
Using formula ()5.16p we calculate the Fourier transformation of the function T (x) 

F[ T (x)] = F[Dh (x)* exp(27rmx)] 

,(m) 



= F[Dh (x)] ■ F[exp(27rmx)] 
(5.19) = r^f^\p) ■ F[^{2TTiax)]. 

By definition of harrow-shaped functions and the Fourier transformation of 6{x — hf3) we have 

F[^{27riax)] = /iexp(27rm/i/3)5(x - /i/3)] 



/i exp(27ri(T/i/3) exp(27rip/i/3) 
/i exp(27ri/i/3((T + p)). 



Hence taking into account (I5.8P we find 

(5.20) (27ri(Tx)] = ^5(o--Fp- /i-^)- 

Taking into account formulas (I5.15|) . (|5.20p from ()5.19p we get 

-1 



F[ T (x)] 



E 



1 



(2vr)2m ^ {p- h-^-ff"^ 



■^6{a + p-h~'P) 



Y,^[a+p-h'^p) 
Y^S{a+p-h-'/3). 



(-1)™ 1 

(2vr)2m 2^ (/1-I/3 - - /i-i7)2™ 



V L 



H -1 



(27r)-™ Y(^"H/3-7)-'t) 



2m 



/3 - 7 / fj/i. 
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From here, since /? and 7 take all integer values, we obtain 

^2m 
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(5.21) 



F[T{x)]=Y,5{a+p-h-^fi) 



(27r)2™Z-(^_^/i)2- 



7 / 



i.e. ah ^ Z. 

Using the identity (jS.Sp . we calculate the inverse Fourier transformation of the function 
Y^5{a +p-h-^P), i.e. 

00 

F-^^5{a +p-h-^P)] = ^ / e~2^"'P(5(a+p-/i-i/3)(ip = ^exp(-2^ix(/i~i/3-fT)) 



(5.22) 



/3 



Now applying the inverse Fourier transformation to both sides of equality ()5.2ip and taking into 
account ()5.22p we have 



T {x) =exp{2Triax) 



1 \m 

^\2m 



(27r)2™ ^ (7 - ahf" 



7 



Hence, going from harrow-shaped functions to functions of discrete argument, we obtain 

1 -1 

-ly 



T[f3] = exp(27riCT/i/3) • 
which completes the proof of Lemma 15. 2[ 



(2vr)2™ Y(7-a/i)^' 



2 m 



ah ^ Z, 



□ 



6. The solution of Problem 3. 
In this section using the results of the previous section we get explicit formula for coefficients 

o 

C([/3]; z) of the lattice optimal interpolation formula and we find unknown constant A. 
Beforehand we give the following result which we use in the proof of the main theorem. 

Lemma 6.1. Let g[/3] be a discrete periodic function, i.e. g[l3\ = g{h(3) = g{hf3 + 7), (3,^ £ Z then 
the following holds 



•1) 



5[/3] = (5[/3]X[o,l][/3])*<^>[/3], 



where X[o 1] [/?] is the discrete characteristic function of the interval [0, 1] and is defined by equality 
Proof. Indeed, using well-known formula (see. [21j ) 
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and periodicity of the function taking into account equality (I5.4p . we have 

k k 

= Y^gm - A;)x[o,i]([/3] -k) = Y,9[P- /i"'^]X[o,i] [/? " h'^k] 

k k 

= Y.Y.9H^mH^[P-i-h-'k] 

k 7 

= 9[l]xm [7] 5^ '^[Z? - 7 - h-^k] = Y 9b]x%i] [7]$[/3 - 7] 
7 fc 7 

= {xmm9m)*m- 

Lemma l6.ll is proved. □ 
The main result of the present paper is the following theorem. 

Theorem 6.2. In the Sobolev space L^2^\o,l) there exists the unique lattice optimal interpolation 
formula of the form U.^) with the error functional coefficients of which have the form 



.2) cimz) = h 



\ 



k 

khiZ 



where L(k) = 



\ -1 

f^2m \ 



I7 — /ifc 



Proof. Applying the operator hD^^\(3]* to both sides of equation (14.30 we obtain 

hD^^\p] * (^B2M * (c([/3],z)x[o,i][/3]) + a) = /il^f^i/?] * B2rn {z - [/?]) , m G (0, 1]. 
Hence taking into account formulas (j5.3p . (|6.ip . (|5.18p we have 

N 

c m,z) -h^c m,z) = hD(^^ m * B^miz - m), m e (o, i]. 

13=1 

By virtue of (|4.4p we find 

Cm,z) = h+ hD^^^ [P] * B2miz -[/?]), [/?] G (0, 1]. 
Hence using the Bernoulli polynomial 

R \RU exp(-27rifc(z - 
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and equalities (I5.17p . (I5.18P we get 



Hence, setting 



/j + /iD;.'"'[/3]* 



exp{—2iTik{z — hf3)) 



{2mk) 



2m 



exp(— 27riA;z) ^(m)rr,^ ., , 



h + {-irh Y: '"^o • (-l)"(2^)^-exp(2vr.fc^/3) 



k 

kh4Z 



(27rfc) 



E 



2?Tt 



(7 - kh) 



2m 



h + hY 

k 

kh0 

h + hY 



(27r)^™exp(27rfA;(/i/3 - z)) 
(27rA;)2™ 

exp(27rik{hf3 — z)) 



^2. 



E 



E 

. 7 



(7 - kh) 

1 



2m 



(7 - kh) 



2m 



.3) 



E 



2m 



(7 - kh) 



2m 



we get (16. 2p . Theorem 16.21 is proved. 



Now using Theorem 16.21 we find A. First, putting the expression (j6.2p of the coefficients (7 
to the left side of equality (14. ip we calculate the following sum 

N 

Q{Z) = Y.Cih]^^)B2m[f3-l] = 
7=1 

/ 



N 



7=1 



Ex - exp(27riA;(/i7 — z)L{k) \ I >. ^ exp(— 27ria(/i/3 — /17) 
+ p;;^ (27ria)2™ 

khfZ / ^ 



V 



Eexp(— 27ria/i/3) 
a^O ^ ' 7=1 

+ E E "lg::'" -E"--pp"^^^ + ^""M. 

A: a^O ^ 7=1 

kh(^Z 
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It is known that 

N 



1 Gxp ( 2vrzo^/i J 



7=1 



(6.4) 



0, ah ^ Z, 

1, ah £ Z, 



N 



E,, , hexT)(2TTi(k + a)h) ■ (1 — exp(2TTi(k + a))) 



7=1 



0, {k + a)h ^ Z, 

1, {k + a)h £ Z. 



(6.5) 

By virtue of equalities (16. 4p . ()6.5I) the expression Q{z) takes the form 



1 -^(^) exp(— 27ri(a/i/3 + /cz) 



when (A; + a)h € Z. 
From here setting 

we have 



{k + a)h = t, t G Z, k + a = th ^, a = th ^ — k, 



_ 1 L(A;) exp(-27rz((t/i-i - A:)/i/3 + kz) 

~ ^ (27rn)2™ + {2mkY^ ^ (t/i^i - A;)2™ 

_ ^ 1 ^ L(A;) ^-^ exp(-27r2((t/3 - + kz) 

1^ (^~)2^+ 2^ (27ri/c)2™ 2^ (t/i-i - /c)2"^ 

El exp(— 27riA;(z — /i/3) • L{k) sr-^ exp(— 27rit/3) 

(27ri7)2'» 2^ (27riA;)2™ 2^ (t/i"! - A;)2™' 

Hence keeping in mind that exp(— 27rzt/3) = 1 and L(A;) = ^ (^fe-i^fc)ijm we obtain 

1 exp(-27rifc(z - hji) 

^ ' > ~ ^ (271^7)2™ + 2^ (2^ife)2- 

Adding and subtracting to the right hand side of equahty ()6.6p the fohowing series 

Eexp(— 27riA;(z — hj3)) 
(27riA;)2™ ' 
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we have 



exp(27rifc(/i/3 - (,_hR\ 



(6.7) 



E 



1 — exp(— 27riA;z) 



Putting the expression (16.7|1 of (^(a:) to the left hand side of equation (I4.ip for A we obtain the following 
expression 



(6i 



exp(— 27riA:z) — 1 



kh€l 



{2mk) 



2m 



Hence clear that when z = hj3 



A = 0. 



Now we show that the expression (16. 2p of the coefficients C{[fi],z) satisfy equality (I4.2p . So, we 
have 



\ ^ i u , exp{2mk{hf3 - z)) 
2^cm,z) = 2^ h + h}_^ L{k) 

13=1 13=1 \ k 

3=1 k 13 = 1 

Nh + hY^ exp(-27rfA;z)L(A:) exp(27ri/i:/i)(l - e^'^*'') 



j^2n 



1-e 



2-Kikh 



1. 



Thus Problem 3 and respectively Problem 2 are solved. 



7. The norm of the error functional 
OF lattice optimal interpolation formulas 

In this section using the results of previous sections we calculate the norm of the error functional £ 
of the lattice optimal interpolation formula. 
The following holds 
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Theorem 7.1. Square of the norm of the error functional of the lattice optimal interpolation 

formula of the form in the space L^2^\o, 1) have the following form 



|Lf)(o,i; 



(-ir 



^ (27riA;)2™ ^ 



2exp{—2Trikz) — 1 



kh& 



{2mk) 



2m 



(7.1) 



L{k) -^-^ exp(27rizt/i 



-1^ 



k 

kh0 



{2mkf'^ ^ {th-^ - /c)2™ 



where L{k) is defined by equality (E 



Proof. From equality (j2.7p when xj. = hk for 



we have 



(-ir 



N 



B2m{0)-2Y,C{h],z) B2m{z-h-f) + 



7=1 



Af N 

+ E C'lM, ^) E Cm,z)B2m{h(3 - hj) 

7=1 /3=1 



Hence taking into account ()4.ip and ()4.2p we obtain 



AT 



N 



{-ir B2m{0)-2^C{b],z)B2m(.Z-hj) + ^ C( , z) (B2m(^ " " A) 

7=1 13=1 
N 

i-ir B2m{0) - C{b],z)B2m{z - h^) - X 
7=1 
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From here using the expression ()6.8p of A we obtain 
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(7.2) 



o 2 



i-iy 



N 



B2m{0)-Y,h[l+ Yl 

7=1 V 



k 

hk6 



exp(27riA;(/i7 — z) L{k) 



exp{—2Tria{z — /17)) exp{—2TTizk) — 1 



(-1) 



2m 



(27rza)2™ 
1 



E 



{2Trik) 



2m 



N 



exp{—2Trikz) L{k) s—^ exp(— 27riaz) 



N 



/ J A;2m 

k a^O 
khiZ 

exp{—2iTikz) — 1 



(27rfa)2™ 



h exp(27rz(a + k)h^) 

7=1 



E 

kh&Z 



{2TTik) 



2m 



From ()7.2p taking into account equahties ()6.4p and (j6.5p we obtain 



(-1)' 



E 



1 



Eexp(— 27riA;z) 



L(A;) 



^2m / > 
khi^Z {k+a)hGZ 



exp{—27riz{k + a)) exp(— 27riA;z) — 1 



(27ria) 



2m 



kh&Z 



{2mk) 



2m 



Whence, setting (A; + a)h = t, a = th ^ — k, we get (j7.ip . Theorem 17. II is proved. 



□ 



Now we show that the expression (|7.ip is zero at the nodes = of the lattice optimal interpo- 
lation formula. 
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= /i/3, then from ()7.ip we obtain 



(-1) 



1 



E 



2exp{-2TTikh/3) - 1 



L{k) v-^ exp(27rit/3) 



kheZ 



{2'Kik) 



2m 



^ (27ri/c)2"^ (27riA;)2™ ^ (27ri/c)2™ • L(/c) 



(27ri/c) 



2m 



(27ri/c)2™ (27riA:)2'» 



0. 



This means that the condition of interpolation is fulfilled and this confirms our theoretical results. 



8. Connection between lattice optimal interpolation formula 

AND optimal quadrature FORMULA IN THE SPACE L^^\{), 1) 

In previous sections we constructed the lattice optimal interpolation formula which has the following 
form 



N 



•1) 



13=1 



where C ([/?], x) are defined by expression (j6.2p . h = jj, N = 2,3, .. 
Integrating equality (|8.ip from to 1, we get 



ip{x)dx 



1 N 



N 



In 



1.2) 



13=1 



\ 



/ h^{hp)dx + h / 



fe 



ex.p{2TTik{h/3 - x)) ■ L{k) 

fc2™ 



The second integral in the right hand side of ()8.2p is equal to zero, since / exp(27rifc(/i/3 — x)dx = 0. 



Then from (j8.2p we get the following well-known quadrature formula 



N 
/3=1 

This is the rectangular formula and optimality of this quadrature formula in the space l'^\o, 1) is 
known (see, for example \23\ 121]). 
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Thus by integrating the optimal interpolation formula in the space L2 (0, 1) we have obtained the 
optimal quadrature formula in the same space. 



9. Numerical results 

First of all here we give explicit formulas for coefficients of lattice optimal interpolation formula 
()1.2p which are very useful in practice. From Theorem 16.21 we get the following 



Corallary 9.1. In the Soholev space (0,1) there exists the unique lattice optimal interpolation 

irror functional coeffi 

2cos{2TT{Nt + j){hl3 - z)) 



formula of the form \1.S\) with the error functional coefficients of which have the form 

\ -In 



(9.1) cm;z) = h 



JV-l 
2 



1 + 

j=i t 



{Nt + jf''^ \ 
when N is odd number and N > m and 

2 cos{27r{Nt + j){h^ - z)) 



E 



1 



cm;z) 



h 



2 ^ 



-1 



i + EE 

j=i t 



(9.2) 



{Nt + j)2'^ 
cos(27r(iVt + iV/2)(/i/3-z)) 



E 



2m 



{Nt + iV/2) 



2m 



E 



1 



(iV7 + N/2) 



2m 



/3 = l,iV 



when N is even number and N > m. 



Proof. Now we simplify the expression ()6.2p of the optimal coefficients C'([/3]; z) 
We denote 



(£>3) 



s = E 



k 



exp(27ri/c(/i/3 — z)) 

J^2m 



E 



1 



-1 



|A^7 - /c|2™ 



Consider two cases. 

The case 1. Let N be odd number, i.e. N = 2M + 1 then k ^ Nt, t e Z. This means 



k 



Nt + 1, Nt + 2,...,Nt + M, 
Nt-1, Nt-2,...,Nt - M. 



Therefore from ()9.3p we have 

M 



s 



EE 

t j=i 



exp{2Tri{Nt + j){h(3 - z)) 



(Nt + j)2m 
exp(27ri(7Vt - j)(/i/3 - z)) 



E 



{N-f - {Nt + j))'^^ 



{Nt-j) 



2m 



E 



(N',-(Nt-j)) 



2m 
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Hence, taking into account that 7 and t take all integer value, replacing 7 by —7 and t by —t in the 
second term of the square brackets we get 

Mr _ , , n / \ -1 

s = EE 

t j=i 

M 



exp{2TTi{Nt + j){hl3 - z)) exp{-2Tri{Nt + j){h(3 - z)) 



EE 

j=i t 



{Nt + j)2m 
2cos(27r(iVt + j)(/i/3-z)) 



{Nt + j)^"^ 
\ -1 



E 



1 



2m 



{Nt + j 



\2m 



E 



1 



(iV7 + 



Taking into account that M = (N — l)/2 and putting the last equality for S into (j6.2p we get (j9.1 
r/ie case 2. Let be even number, i.e. = 2M then k 7^ A^i, t G Z. This means 



iVt + 1, Nt + 2, iVt + M - 1, TVi + M, 
iVt - 1, Nt- 2, iVt - M + 1. 



Therefore from ()9.3p we have 

M-l 

s = EE 



exp(27ri(iVt + j)(/i/3-z)) 



{Nt + j 



\2m 



E 



+ 



exp(27ri(iVt - - z)) 



(Ar^_(iVt + j)) 



2m 



-1 



E 



(A^7 - {Nt - j))^™ 



-E 



exp(27ri(A^t + M){hp - z)) ^ exp{2m{Nt - M){hj3 - z)) 



{Nt + M) 



2m 



{Nt - My^ 



E 



{N-f + 7^)2"^ 



Hence as in the case 1, taking into account that 7 and t take all integer value, replacing 7 by —7 and 
t by — t in the second terms of the square brackets we get 



S 



M-l 

EE 

t j=i 

t 

M-l 



exp{27ri{Nt + j){h/3 - z)) ^ exp{-2Tri{Nt + j){h/3 - z)) 



{Nt + j 



\2m 



{Nt + j 



\2m 



E 



1 



-1 



(A^7 + j)^™ 



+ 



EE 

j=i t 



exp{2Tri{Nt + M){h/3 - z)) ^ exp{-27ri{Nt + M){hp - z)) 
{Nt + M)2™ 

2cos{2Tr{Nt + j){hl3 - z)) 



+E 



{Nt + j)2'» 
cos{2Tr{Nt + M){h/3 - z)) 



E 



1 



{Nt + M)2™ 
\ -1 



E 



(A^7 + M)2™ 



{Nt + M) 



2m 



E 



(A^7 + j)2m 



1 

2m 



{Nj + M) 



Taking into account that M = N/2 and putting the last equality for S into (j6.2p we get (j9.2 
Corollary 19.11 is proved. 



□ 
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Here in numerical examples we consider interpolation of two 1-periodic functions: 

1) (pi{x) = sin(27rx), 

2) the Bernoulli polynomial of degree 10, i.e. 

M^) = B,o{x) = - + 15^8 _ + _ 3 ^ ^ 5 

2 2 bb 

We interpolate these two functions with optimal interpolation formula 

5 

(9.4) P^(x) = ^(7([/3];x)(^(/i/3) 

(3=1 

for m = 1,2,3,4. For simplicity we have taken N = 5. Since in our case N = 5, h = 1/5 as optimal 
coefficients we use formula (|9.ip . One can use formula (j9.2p for even number N of the nodes. 

Below in each case graphs of the optimal coefficients and graphs of absolute errors between optimal 
interpolation formula (j9.4|) and functions sin(27rx), Biq{x) are respectively given. 

The case m = 1. 




Figure 1. Graphs of optimal coefficients for m = 1, = 5: a) C([l],x), b) C([2],x), 
c) C([3],x),d) C([4],x),e) C([5],x). 




Figure 3. Graphs of optimal coefficients for m = 2, iV = 5: a) C([l],x), b) C([2],x), 
c) C([3],x),d) C([4],x),e) C([5],x). 
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Figure 4. Graphs of absolute errors for m = 2, A'' = 5: a) | sin(27rx) — P^{x)\, b) 
\Bioix) - P^{x)\. 



The case m = 3. 




Figure 5. Graphs of optmial coefficients for m = 3, A?" = 5: a) C([l],x), b) C([2],x), 
c) Ci[3],x), d) Cm,x), e) Ci[5],x). 




Figure 7. Graphs of optimal coefficients for m = 4, A?^ = 5: a) C([l], x), b) C{[2],x), 
c) C([3],x), d) C{[4],x), e) C{[5],x). 
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Figure 8. Graphs of absolute errors for m = 4, = 5: a) | sin(27ra;) — P^{x)\, b) 
\Bio{x) - P^{x)\. 

From Figures [2|4|6|8I we can conclude that the absolute errors between given functions and optimal 
interpolation formula is decreasing as m is increasing. 
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